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Abstract. We report high resolution (l"-2") imaging of the CO 2-1 line and the millimeter continuum in the 
proto-planetary nebula He 3-1475. The observations reveal the presence of a massive (~ 0.6 Mq) envelope of 
molecular gas around the origin of the remarkable bipolar jet system seen in optical images with the HST. The 
CO kinematics are well modeled by an expanding, bi-conical envelope: the prominent, high- velocity ('^ 50 kms~^) 
wings seen in single-dish CO spectra arise where the sides of the bi-cones are projected along the line of sight. The 
continuum is detected at 1.3 mm and 2.6 mm and is due to thermal emission from warm ('^ 80 K) circumstellar 
dust. The structure, kinematics, and expansion time of the envelope provide strong evidence for entrainment of 
the molecular gas by the high velocity jets. The observations support an evolutionary scenario in which a period 
of enhanced mass loss by the central star is followed by the development of the bipolar jets which burst through 
the molecular envelope. The jet-envelope interactions play a crucial role in shaping the subsequent ionized nebula. 
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1. Introduction 

Collimated, bipolar outflows or "jets" are an important 
feature of the early evolution of planetary nebulae (PNe) 
that has only recently been widely recognized. Many 
proto-PNe and PNe are now known to show evidence of 
bi-polar or multi-polar structures, or point symmetries, 
which have been produced by the action of symmetric 
jets from the central star (see, e.g., Kastner et al. t20g0|, 
and these structures are sufficiently common that possibly 
most or all PNe pass through this phase (Sahai & Trauger 

The jets are most active in the early phases of PNe 
formation, and their effects can be detected in proto-PNe 
as high velocity wings in low angular resolution, molecu- 
lar line spectra of the neutral circumstellar envelopes (e.g., 
Cernicharo et al. 1989, Young et al. 1992; Bujarrabal et al. 
I2001|l . In the few cases observed at high angular resolution, 
the wings are seen to arise in directed outflows of entrained 
molecular gas (e.g., Cox et al. 2000, 2 003 , 1 . The origin of 
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the jets that cause the outflows is not well understood, but 
the jets clearly have major effects on the structure and dy- 
namics of the neutral circumstellar envelopes from which 
the ionized nebulae form (Huggins et al. 1996), and conse- 
quently they play a key role in the early shaping of PNe. 
Well studied examples of young PNe that show evidence 
for interactions between the jets and the neutral circum- 
stellar gas include BD-h30°3639 (Bachiller et al. IMlUjl . 
KjPn 8 (Forveille et al. n^Sjl and Ml-16 (Huggins et al. 

He 3-1475 (IRAS 17423-1755) is an extreme exam- 
ple of a proto-PN with highly collimated bipolar outflows, 
first discussed by Riera et al. H1995|l and Bobrowski et 
al. I|1995|l . It is of special interest because of the high ve- 
locities and remarkable structure seen in high resolution 
optical observations made with the HST (Borkowski et al. 
'1997": Borkowski & Harrington l2001l Sanchez-Contreras & 
Sahai 200lj Riera et al. 12003(1 . The outflows terminate in 
a series of knots that are point symmetric about the cen- 
tral star, and indicate formation by episodic jets whose 
direction is time dependent. 

Little is known about the neutral circumstellar matter 
in He 3-1475. IRAS fluxes indicate the presence of a cir- 
cumstellar dust envelope (Parthasarathy & Pottasch ll989|l 
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Fig. 1. CO spectra in the 1-0 (115 GHz) and 2-1 (230 GHz) lines towards He 3-1475, obtained with the IRAM 30 m 
telescope. 
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Fig. 2. Millimeter maps of He 3-1475 obtained with the IRAM interferometer. Left: the continuum at 2.6 mm. Center: 
the continuum at 1.3 mm. Right: the velocity integrated CO 2-1 line. The contour intervals are 1 mJy beam~^ (2.6 mm), 
2 mJy beam~^ (1.3 mm), and 50 K kms~^(CO); dashed contours are negative. The field center is 17'M5™14!17, 
— 17°56'47'.'0 (J2000), and the beam size is shown in the lower right of each panel. 



which is also seen in optical images as a dark lane cross- 
ing the nebular axis to the south-east of the central star. 
OH maser emission at 1667 MHz has been observed by 
te Lintel Hekkert H1991|l and Bobrowski et al. (1995), and 
CO spectra have been reported by Knapp et al. ( 1995 ) and 
Bujarrabal et al. (2001). In this paper we report high an- 
gular resolution observations of the CO emission and the 
millimeter continuum to study the relation of the neutral 
circumstellar gas to the outflows. 

2. Observations 

High sensitivity, single-dish observations of He 3-1475 were 
made in the 2.6 mm CO J = 1 — (115 GHz) and 1.3 mm 
J = 2 - 1 (230 GHz) fines using the IRAM 30 m tele- 



scope at Pico Veleta, Spain. The observations were made 
in September 1997, using 512 x 1 MHz filterbanks. The 
half power beam size of the telescope is 11" and 22" at 
the frequencies of the 2-1 and 1-0 lines, respectively. The 
calibration was made using tfie chopper wheel technique, 
and the line intensities are reported fiere as main beam 
brightness temperatures. 

High angular resolution observations were made in the 
CO 1-0 and 2-1 lines and in the nearby continuum dur- 
ing February and March 1998 using the IRAM interfer- 
ometer at Plateau de Bure, France. The array consisted 
of five 15 m antennas, equipped with SIS heterodyne re- 
ceivers. The observations, centered on He 3-1475, were 
made with two configurations of the array, with baselines 
up to 280 m. The primary beam size of the interferometer 
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Table 1. CO observations of He 3-1475 made with the IRAM 30 m telescope 
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Fig. 3. CO 2-1 map of He 3-1475 superposed on the HST 
WFPC2 image in [N ii] A6584 (Borkowski et al. lTMTji . The 
CO contours are the same as in Fig. 2. 



is 22" at 1.3 mm and 44" at 2.6 mm. The effective veloc- 
ity resolution of the line observations used for the analy- 
sis is 8 kms""'^. The continuum observations were made at 
frequencies of 115.3 and 231.5 GHz, with effective band- 
widths of 320 and 640 MHz, respectively. 

The RF passband and amplitude were calibrated us- 
ing 3C273 and the phase calibration was performed ev- 
ery 20 minutes using J1730-130 and J1830-210. The 
uv data were Fourier transformed and CLEANed, us- 
ing the Clark algorithm and the restored Gaussian clean 
beam is 3"4 x 1"6 (PA = 177°) at 2.6 mm for both 
the line and continuum, and 1"5 x 0"7 (PA = 8°) and 
2"2 X 0"8 (PA = 0°) at 1.3 mm for the hne and contin- 
uum, respectively. The adopted field center for the maps 
is 17^45'"14!17, -17°56'47"0 (J2000). 



3. Results 

The CO 1-0 and 2-1 spectra obtained with the 30 m tele- 
scope towards He 3-1475 are shown in Fig. 1. A small map 
was made with 5" spacing around the center position, but 
the molecular emission was not found to be extended with 
respect to the telescope beam. The parameters of the lines, 
based on Gaussian fits to the spectra, and the correspond- 
ing line fluxes are given in Table 1. 



The CO 1-0 and 2-1 lines were both detected at high 
angular resolution with the interferometer, and the obser- 
vations recover 97% and 72% of the single-dish fluxes for 
the 1-0 and 2-1 lines, respectively. The observed distri- 
bution of CO is shown in the velocity integrated map of 
the 2-1 emission in the right hand panel of Fig. 2. The 
position of the peak CO emission is at ±0"15 ± 0"15, 
±0'.'05 ± 0'.'15 relative to the map center, and the emis- 
sion is extended with respect to the telescope beam, with 
a deconvolved source size of - 1"6 (FWHM). The CO 
distribution is compared with the optical structure of the 
nebula in Fig. 3, and the CO kinematics are shown in 
channel maps and velocity-position maps, together with 
model simulations, in Figs. 4-6. The interferometer obser- 
vations in the CO 1-0 line produced essentially the same 
results as the 2-1 observations, but with a factor of two 
lower resolution (because of the longer wavelength) , so are 
not discussed further. 

The millimeter continuum of He 3-1475 was detected at 
1.3 mm and 2.6 mm with the interferometer, and maps of 
the emission are shown in the center and left hand panels 
of Fig. 2. At 2.6 mm the emission is not extended with 
respect to the telescope beam, and the measured flux is 
5.3 ± 1.2 mJy. At 1.3 mm the observed emission (2" 8 ± 
0"2x l'.'l±0"4, PA = 12°) is slightly extended with respect 
to the beam, and the measured flux is 31 ± 4 mJy. The 
position of the peak emission in the 1.3 mm continuum is 
at +0'.'04 ± 0"05, +0"40±0"12 relative to the map center, 
and is essentially coincident with the position of peak CO 
emission, within the uncertainties. 

4. Properties of the envelope 

4.1. Overview 

The CO spectra of He 3-1475 (Fig. 1) are much broader 
than the 5-15 kms^^ lincwidths typically seen in AGB en- 
velopes, and the shapes of the profiles are different. They 
show prominent, high velocity wings with no clear low 
velocity component corresponding to an undisturbed en- 
velope: thus most or all of the molecular gas appears to 
participate in the high velocity flow. 

From the measurements given in Table 1, the systemic 
LSR radial velocity (Vo) of the molecular gas is well deter- 
mined to be T4 = 47.5 ± 0.8 kms^^. This corresponds to 
a heliocentic velocity of 34.5 kms^^, and is in good agree- 
ment with the (heliocentric) velocity of 37.6 ± 1.2 kms~^ 
recently reported by Borkowski & Harrington H2001|l . 
based on five stellar lines. 
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Fig. 4. Channel maps of the CO 2-1 emission observed in He 3-1475. The channels are 8 kms~^ wide, and are centered 
at the velocities given in the upper left of each panel. The contour interval is 0.5 K; the dashed contours are negative. 
The arrows indicate the major (jet) axis (PA = 135°), and the beam size is shown in the lower right panel. 




Fig. 5. Theoretical channel maps of the CO 2-1 emission for the biconical model discussed in the text. The channels 
are 8 kms~^ wide, and are centered at the velocities given in the upper left of each panel. The contours are from 5% 
to 95% (in steps of 10%) of the peak emission. 
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Fig. 6. Position- velocity maps of the CO 2-1 emission in He 3-1475 along the major {PA ~ 135°) and minor {PA = 
45°) axes. Left: observed data; the contour interval is 0.5 K; dashed contours are negative. Right: theoretical maps for 
the biconical model discussed in the text. The contours are from 5% to 95% (in steps of 10%) of the peak emission. 



The comparison of the CO map and the HST WFPC2 
image in [N ii] A6584 (Borkowski et al. [TM7| in Fig. 3 
shows that the CO emission peaks close to the center of 
the nebula, but is quite limited in extent compared to the 
optical image. The molecular emission covers the central 
dark lane and extends out along the opening of the bipolar 
structure, but does not envelop the extended, coUimated 
jet system. The characteristic radius of the CO map (O'.'S) 
corresponds to 7 x 10^^ cm at a distance of 5.8 kpc, which 
we adopt for the distance to He 3-1475 (Riera et al. 1^003,) . 

4.2. CO kinematic structure 

Although the CO emission is not very extended with re- 
spect to the telescope beam, the velocity-resolved CO 
observations provide insights into the overall kinematic 
structure of the molecular gas. The complete data cube is 
shown in Fig. 4 in the form of channel maps. These exhibit 
an approximate position-velocity symmetry about the 
central position and velocity, e.g., the emission in the ex- 
treme blue shifted channels (at —6, +2, and +10 kms^^) 
is offset to the NW of the field center and that of the ex- 
treme red channels (-1-90, -1-98 and -1-106 kms~^) is offset 



to the SE. The maps also exhibit a rough axial symmetry 
about the major axis of the optical nebula (the jet axis) 
at PA = 135°, although the CO intensities are affected by 
the shape of telescope beam, which is elliptical and lies 
at an angle of 53° with respect to the nebula axis. Thus 
the arc-like features in blue channels -1-10 and -1-18 kms""'^ 
and red channels -1-82 and -1-90 kms^^, which point away 
from the center, appear stronger where the emission lies 
along the primary axis of the beam (see §4.3 below). 

The velocity-position maps through the data cube 
along the major and minor axes shown in Fig. 6, demon- 
strate that the whole CO structure is tilted toward us 
(blue shifted) to the NW and away from us (red shifted) 
to the SE along the major axis. This tilt is in the same 
sense as that of the optical bipolarity and jets (e.g., Riera 
et al. 2003). In addition, the distribution of intensity in the 
map along the major axis indicates that the molecular gas 
forms an expanding, open ended, bi-conical structure that 
can also be traced in the channel maps in Fig. 4. Thus 
the extreme, blue shifted emission (in channels —6 and 
+2 kms^^) arises in the lower rim of the cone facing to- 
ward us; this rim lies nearly in the direct line of sight, and 
gives rise to the highest blue shifted gas seen in the wings 
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of the single dish spectra. At intermediate, blue shifted 
channels (10-26 kms~^), the arc-like structure becomes 
more prominent as the channel includes a larger cross sec- 
tion of the cone, and closer to the systemic velocity (at 
-f42 kms^^) the channel also includes the upper rim of 
the cone. At red shifted channels, the overall symmetry is 
reversed for the cone facing away from us. 

4.3. Model of the CO kinematics 

In order to quantify our interpretation of the CO data, 
we have constructed a simple, biconical model of the CO 
emission for comparison with the observations. The model 
is characterized by the radius at the equator where the 
cones intersect {Ri), the radius of the open ends (i?o), the 
height of each cone (h), and the inclination of the sym- 
metry axis to the line of sight (i). We assume a constant 
thickness {t) , and a relative density within the cones that 
varies as a power law from the center (p ~ r"). The emis- 
sion is taken to be optically thin (see below), and the 
expansion is assumed to be homologous, with the radial 
velocity v ^ r. This velocity law is motivated by the ob- 
served structure in the data cube, and by similar ballistic 
flows seen in other proto-PNe (e.g., AFGL 618, Cox et al. 

mm- 

From comparison with the observations, the parame- 
ters of a best-fit model are Ri = 0'.'45, Ro = 0'.'92, h = 1'.'5, 
t = 0'.'5, n — —0.5, i — 40°, and a velocity gradient of 
31.4 kms^^ arcsec^^. The linear dimensions are given in 
arc seconds on the sky, since they scale with distance. 
The thickness of the walls of the bicones is not actually 
resolved by the observations, and is set to a nominal value. 
The parameter n is constrained by the relative intensity of 
the extended CO emission to that near the center, and the 
dimensions of the bicones, the velocity gradient, and the 
inclination are jointly constrained by the map sizes at dif- 
ferent velocities, and the relative projections of the sides 
of the cones in the position-velocity map. For example, a 
general constraint on the inclination and the opening an- 
gle of the bicones is provided by the major axis position- 
velocity map, where the blue shifted, lower rim of the cone 
facing toward us is very close to the line of sight, and the 
upper rim is also blue shifted. In this case, the inclination 
angle of the symmetry axis i ^ 45° and the opening angle 
of the cone from the center is ^ 45° (otherwise the upper 
rim would be red shifted). Similarly the inclination cannot 
be much less than ~ 40° or the upper rim would project 
to higher velocities than are observed. 

The results of the best-fit model are shown in Fig. 5 
and Fig. 6 (right panel). In spite of the extreme simplic- 
ity of the model, it can account for the main features of 
the observations. The inclination of the CO bicones to the 
line of sight is essentially the same as that of the optical 
jets (Borkowski & IIarrine:ton ,2001 ). and the biconical ge- 
ometry indicates that molecular gas lies around the base 
of the jet flows. As noted above we do not resolve the 
walls of the bicones, but the parameter n together with 
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Fig. 7. The infrared-radio spectrum of IIe3-1475. The 
flux measurements are from IRAS (squares), this paper 
(crosses), and Knapp et al. I|1995|l (circle). The smooth 
curve is for thermal emission from dust at a temperature 
of 81 K. 

the assumed constant wall thickness determines the mass 
distribution in the bicones: for the solution n = —0.5, the 
mass per unit length projected along the major axis is ap- 
proximately constant, which is reasonable for such a flow. 

The kinematic model is axi-symmetric but our calcula- 
tions of the simulated observations also include the ellip- 
tical telescope beam which is not aligned along this axis. 
The effects can be seen in Figs. 5 and 6. They include 
the asymmetries in intensity of the arc-like features in the 
channel maps noted earlier, and a more subtle effect in 
a shift from red to blue (from NE to SW) in the minor 
axis position-velocity map, which is also apparent in the 
real data. Two features of the observations are not re- 
produced by the model simulations. The first is a slight 
channel-to-channel shift in the position of the emission at 
high velocities, and this could be accommodated by some 
degree of curvature in the sides of the bicones. The sec- 
ond is slightly higher equatorial emission in the observa- 
tions than in the model, which could be accommodated by 
higher densities or a modified geometry near the equator. 
However, the emission is probably at least partially thick, 
especially near the equator, and we have not included this 
in the simple model, so fine tuning with additional param- 
eters would be unwarranted. Nevertheless, independent of 
optical depth effects, the basic kinematic structure is well 
accounted for by the simple model. 

4.4. The circumstellar dust component 

Our millimeter continuum observations fill a large gap in 
the observed long wavelength spectrum of He 3-1475, and 
they provide an estimate of the mass of circumstellar dust. 
As shown in Fig. 7, the only other long wavelength obser- 
vations of He 3-1475 arc from IRAS, and a single radio 
observation at 3.6 cm. The centimeter continuum is likely 
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to be free-free emission from the compact ionized core, 
and from comparison with the flux in H a Borowski et al. 
H1995|l suggest that the emission is optically thick. An op- 
tically thick radio spectrum which varies as does not, 
however, extend into the millimeter region because the 
115 GHz (2.6 mm) flux would then be much higher than we 
observe. The free-free emission spectrum likely becomes 
optically thin and levels off between 10 and 100 GHz. 

The millimeter flux increases quite steeply with in- 
creasing frequency, and from the overall shape of the spec- 
trum the main contribution to the observed 1.3 mm and 
2.6 mm flux is probably the long wavelength tail of the 
dust emission seen in the infrared. This is supported by 
the fact that we are able to simultaneously provide a good 
fit to the 25-100 ^m (color corrected) IRAS fluxes and to 
the observed millimeter fluxes, with a single temperature 
dust component. Formal flts, varying both the tempera- 
ture and the index of dust emissivity p (where Q{iy) ^ v^), 
give Td = 81 ± 4 K, and p = 0.99 ± 0.07. 

With this value for the temperature, we estimate the 
mass of the cool dust component using the optically thin 
expression — F^D^ /{xuBi^(T)), where is the flux, 
Xu the emissivity per unit mass, and B,^ the Planck func- 
tion at frequency v, and D is the distance. Using the 
IRAS flux at 60 ^im, xeo = 150 cm^ g-^ (Jura 1986), 
and D = 5.8 kpc, we find Md = 6.4 x 10"^ Mq. A second 
component with a higher temperature and a much smaller 
mass, which we ignore, can account for the additional, 
short wavelength emission in the spectrum in Fig. 7. 

4.5. Mass of circumstellar gas 

The mass of circumstellar gas can be estimated from the 
dust mass given above, by adopting a value for the gas-to- 
dust ratio. For transition objects with high mass loss rates, 
the gas-to-dust ratio appears to be somewhat smaller than 
for typical AGB stars, and we adopt the value of ^ 100 
from Knapp et al. H1993|l . This yields a mass of circum- 
stellar gas Mg = 0.64 Mq. 

The CO observations provide a second estimate for the 
mass of circumstellar gas. In this case, it is a lower limit 
because the low (- 0.9) CO 2-1/1-0 flux ratio (Table 1) 
suggests that the CO lines may be at least partially op- 
tically thick, although an alternative possibility is that 
the CO is sub-thermally excited. In any event we ob- 
tain a lower limit to the mass of molecular gas using 
the optically thin formula given by Huggins et al. (^5^. 
For a distance of 5.8 kpc and a representative CO/H2 
abundance of 2 x 10^* which is commonly assumed for 
oxygen-rich envelopes (e.g., Kahane & Jura ll99Hl . we find 
Mg <; 0.19 Mq, consistent with the estimate given above. 

These values are in accord with a further estimate of 
the mass of circumstellar gas, given by Bujarrabal et al. 
H2001|l based on observations of the ^^CO lines. This ap- 
proach minimizes the effect of line opacity because the 
lines are likely to be optically thin, but relies on an as- 
sumed value for the ^^C0/H2 abundance. For their as- 



sumed value (2 x 10 and our adopted distance of 
5.8 kpc to He 3-1475, Mg = 0.85 Mq, consistent with 
the above values. 

5. Jet-envelope interactions 

Our observations of the molecular emission in He 3-1475 
support an evolutionary scenario in which a period of en- 
hanced mass loss by the central star is followed by the 
development of bipolar jets that burst through the sur- 
rounding molecular gas. 

5.1. Mass loss rate of the precursor 

The current mass of circumstellar molecular gas around 
He 3-1475 is substantial (§4.5), and it was presumably 
ejected by the star at a moderate velocity characteristic 
of the AGB. For a velocity of ~ 15 kms^^, the size of 
the CO map that we observe (§3) implies an ejection time 
scale of ^ 1,500 yr, and a corresponding mass loss rate 
over this time of ^ 1.3 x lO^'' MQyr~^, using our lower 
limit on the CO mass. This mass loss rate is much larger 
than typically seen on the AGB, but is characteristic of 
some other transition objects (e.g., AFGL 2688, Jura et 
al. I2000|l , and is either an intrinsic part of the final evo- 
lution of single stars on the AGB, or the resiflt of binary 
interactions. 

5.2. Envelope entrainment 

In addition to enhanced mass loss, the kinematics and 
structure of the circumstellar gas indicate that most or 
all of the molecular envelope has been affected by recent 
interaction with the jets. First, the CO velocities (up to 
^ 50 kms^^) are significantly larger than the expansion 
velocities of AGB stars, but are much less than velocities 
seen in the ionized gas close to the jet axes (Borkowski & 
Harrington 120011 Sanchez-Contreras & Sahai I^UITIl Riera 
et al. l2003|l . Second, the observations show that the molec- 
ular gas forms an expanding bi-conical structure around 
the base of the optical bipolar flows. Thus both the struc- 
ture and kinematics provide strong evidence for entrain- 
ment of the molecular gas. 

The kinematic time scales support this view. From 
proper motion studies using HST, the expansion time 
scales {r/v projected on the plane of the sky) of the ion- 
ized knots at ±6" from the center are Texp ^ 450-550 yr, 
and for the most distant knots at ±7'.'5, Toxp ~ 600 yr 
(Riera et al. 12003. Borkowski & Harrington ^UU^. For 
comparison, the kinematic time scale of the molecular gas 
is obtained by combining the (angular) velocity gradient 
of our kinematic model (§4.3) and the adopted distance of 
5.8 pc, which give Tcxp ^ 875 yr. Thus the CO outflows 
and the earliest jets are comparable in age. 

There are large velocity gradients from the axes of the 
jets to the sides of the CO bicones, as expected if the jets 
burst through the envelope. The maximum CO velocity 
that we observe is ~ 50 kms^^, which is in fact close to 
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the maximum velocity for molecular gas to survive accel- 
eration in a single shock (Draine et al. I1983|l . although 
it is possible for molecules to reform at higher velocities 
in post-shocked gas (HoUenbach & McKee [19891 Neufeld 
& Dalgarno 1989). Details of the interaction need further 
study, but it seems likely that much of the fast, ionized 
material seen in the optical outflows has been produced by 
the destruction of molecules in the entrainment process. 
Similarly, there may also be an intermediate component 
of neutral atomic gas produced by the interactions which 
could be detected in species such as C ii, C i, and O i. 

Molecular H2 has been detected in the axial knots 
themselves (Harrington et al. l2(J0(Jfl . This is probably at 
very high velocity (although no spectroscopy has yet been 
done) and the emission is likely formed in dense, post- 
shocked gas, possibly under conditions similar to those in 
AFGL 618 discussed by Cox et al. ijW)^ . 

Weak OH maser emission has also been observed to- 
ward He 3-1475 (te Lintel Hckkcrt lTMTl Bobrowski et al. 
I1995|l and is distributed in multiple maser spots within 
~ O'.'S of the center, with a velocity range ^ 20-75 kms~^, 
roughly centered on the CO systemic velocity. The kine- 
matic structure of the OH masers in Fig. 5 of Bobrowski 
et al. H1995|l appears to show a strong velocity gradient 
across the envelope which has been remarked on by oth- 
ers, e.g., Riera et al. H2()()3|l . This apparent gradient is, 
however, not real. It is an artifact of ascribing positive and 
negative radial offsets to the red and blue components, and 
the actual distribution of the maser components shows no 
clear cut geometry (Zijlstra et al. I2()()l|l . Given the OH 
maser velocities, the spots are probably located in the in- 
ner regions of the CO bicones near the equatorial plane. 

5.3. Envelope dynamics 

It has long been suspected that stellar radiation pressure 
may not power the highest mass loss rates in post-AGB 
stars (e.g., Knapp et al. ll982|l . and recent survey work by 
Bujarrabal et al. (20(nj has quantified this for the molecu- 
lar outflows seen as high velocity wings in single dish CO 
spectra. We follow their approach in estimating the lin- 
ear momentum (P) and energy (E) of the outflow using 
the observed values along the line of sight, and correcting 
for the inclination. We use the position-velocity map of 
Fig. 6 to define the outflow as the emission at velocities 
|F - Fo| > 5 kms"^ For d = 5.8 kpc and i = 40° adopted 
earlier, we find P ^ 1 x lO'^^ gcms~^ and E ^ 2x 10'^^ erg 
using the CO emission. If we use the mass estimate based 
on the millimeter continuum and assume a velocity dis- 
tribution as in the CO profile, the actual values are ^ 3 
times larger than these limits. The results are similar to 
those obtained by Bujarrabal et al. H2001|l for He 3-1475, 
based on the ^^CO lines. 

The luminosity of He 3-1475 at a distance of 5.8 kpc 
is L = 12,600 Lq (Riera et al. I2003|l . If radiation pres- 
sure from the star drives the molecular outflows, the time 
scale to generate the observed momentum (defined by the 



equation t{L/c) = P) is ^ 2 x lO'* yr. However, as dis- 
cussed in §5.2, the CO expansion time is only ~ 875 yr, 
and the actual time to accelerate the molecular gas to the 
observed velocities is probably less. Thus radiation pres- 
sure is unable to drive the flow by a large margin, as found 
by Bujarrabal et al. (I200ip for more than 20 cases. This re- 
sult is not unexpected in He 3-1475 because the geometry 
and kinematics discussed above provide strong evidence 
that the molecular gas is accelerated by entrainment in 
the jets, and these highly coUimated structures are not 
likely generated by radiation pressure. 

In He 3-1475 the jets appear to have easily penetrated 
the dense circumstellar gas close to the star: the molecular 
gas is peripheral to the outflows, and is at relatively low 
velocities compared to the jet axes. In these circumstances, 
it might be expected that the energy and momentum of 
the flow estimated from the molecular gas alone would be a 
relatively small part of the total, i.e., that the values given 
above are strong lower limits. It is therefore surprising that 
the energy and momentum in the fast outflows estimated 
by Riera et al. (j2003|l from optical observations are P ~ 
10"^^ gcms~^ and E ~ 10^^ erg, which are less than the 
values we estimate from the molecular emission, the value 
of P less by a factor of 100. 

This discrepancy requires explanation and we offer two 
solutions. One is that the material close to the jet axes is 
largely neutral. Riera et al. ( 2003} have estimated the flow 
parameters based on a density determined from [S 11] line 
ratios. If the material is largely neutral, as expected from 
calculations of jets in YSOs (e.g., Safier [TM^ Shang et 
al. 12002(1 . then the actual densities could be much larger, 
and the dynamical quantities correspondingly higher. A 
second possibility is that the primordial jets are in fact 
components of a wide angle wind. If so, the wide angle 
wind might couple well to the molecular gas, with only 
the central, highest velocity material breaking through the 
envelope to form the visible jets. These scenarios need to 
be examined in other proto-PNe. 

6. Conclusions 

The millimeter observations reported here provide basic 
information on the properties of the neutral circumstel- 
lar matter around the remarkable proto-planetary nebula 
He 3-1475. The observations also underscore the impor- 
tance of neutral circumstellar gas in the early development 
of PNe. 

He 3-1475 is surrounded by a massive circumstellar 
envelope (kI 0.2 Mq and ^ 0.6 A/q from our CO and 
continuum observations, respectively), which has only re- 
cently been ejected by the central star at a high mass loss 
rate (^ 1 x 10^** Mq yr^^). The structure and kinematics 
of the CO emission are well modeled with an expanding 
bi-conical envelope, and lead to the conclusion that the 
molecular gas has been entrained in the sides of the jets 
as they burst through the envelope. The expansion time 
scales of the CO emission and the jets support this view. 
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Although He 3-1475 is an extreme object on account 
of the very high velocity of the jets and its well developed 
point symmetric structure, the evolutionary scenario out- 
lined above is very similar to other newly forming PNe we 
have studied at high resolution. AFGL 2688 (Cox et al. 
innnOJ, CRL eiS (Cox et alHrOS); Ml-16 (Huggins et al. 
EOnnj, KjPn 8 (Forveille et al. 11998.1 . and NGC 7027 (Cox 
et al. I2002|l form an approximate evolutionary sequence 
in which the ionized nebula turns on and becomes more 
dominant, and in each case there are prominent multiple 
jets, or single jets which have changed direction. He 3- 
1475 clearly belongs with this class. Besides constraining 
the origin of the jets and the physics of jet-envelope inter- 
actions, the observations of these objects demonstrate the 
importance of jets in the early shaping of the neutral cir- 
cumstellar envelopes which play a key role in determining 
the morphology of the mature PNe. 
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